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Abstract New measurements using radio and plasma-wave instruments in in- 
terplanetary space have shown that nanometer-scale dust, or nanodust, is a 
significant contributor to the total mass in interplanetary space. Better measure- 
ments of nanodust will allow us to determine where it comes from and the extent 
to which it interacts with the solar wind. When one of these nanodust grains 
impacts a spacecraft, it creates an expanding plasma cloud, which perturbs the 
photoelectron currents. This leads to a voltage pulse between the spacecraft body 
and the antenna. Nanodust has a high charge/mass ratio, and therefore can be 
accelerated by the interplanetary magnetic field to speeds up to the speed of the 
solar wind: significantly faster than the Keplerian orbital speeds of heavier dust. 
The amplitude of the signal induced by a dust grain grows much more strongly 
with speed than with mass of the dust particle. As a result, nanodust can produce 
a strong signal, despite their low mass. The WAVES instruments on the twin 
Solar TErrestrial RElations Observatory spacecraft have observed interplanetary 
nanodust particles since shortly after their launch in 2006. After describing a new 
and improved analysis of the last five years of STEREO /WAVES Low Frequency 
Receiver data, a statistical survey of the nanodust characteristics, namely the rise 
time of the pulse voltage and the flux of nanodust, is presented. Agreement with 
previous measurements and interplanetary dust models is shown. The temporal 
variations of the nanodust flux are also discussed. 
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1. Introduction 

Radio and plasma-wave instruments have measured interplanetary dust in situ 
since 1983. Dust grains were first detected with electric sensors in planetary 
environments, during the G-ring crossing of Saturn by Voyager 2 (Aubier, Meyer- 
Vernet, and Pedersen, 1983; Gurnett et al., 1983). Dusts have been measured 
in situ by such instruments in planetary environments, during the ring-plane 
crossings of Uranus (Meyer- Vernet, Aubier, and Pedersen, 1986; Gurnett et al., 
1987) and Neptune (Pedersen et al, 1991; Gurnett et al, 1991), the E-ring 
crossing of Saturn by Voyager 1 (Meyer- Vernet, Lecacheux, and Pedersen, 1996), 
and recently by Cassini (Kurth et al, 2006). Dusts have also been observed in 
comctary environments, with the International Cometary Explorer spacecraft, 
which crossed the tail of comet Giacobini-Zinner in 1985 (Gurnett et at, 1986). 
Later, the Vega spacecraft measured dust particles near comet Halley (Oberc, 
Parzydlo, and Vaisberg, 1990; Oberc and Parzydlo, 1992; Oberc, 1993). 

In addition to the measurements near outer planets by classical cosmic dust 
analyzers (see Hsu, Kriiger, and Postberg (2012) and references therein), nano- 
sized dust grains have been measured recently in the vicinity of Jupiter by 
Cassini (Meyer- Vernet et al., 2009a), and in the solar wind at 1 Astronomical 
Unit (AU) (Meyer- Vernet et al, 2009b; Zaslavsky et al, 2012; Meyer- Vernet and 
Zaslavsky, 2012). 

In this article, we present an improved analysis of the detection of nanome- 
ter dust grains by the Solar TErrestrial RElations Observatory /WAVES Low 
Frequency Receiver (LFR) (Section 2), from whose measurements we deduce 
nanodust characteristics (Section 3). We study five years of data between 2007 
and 2011. This data set corresponds to 3 980 919 spectra for STEREO-Ahead 
(A) and 3 952 952 for STEREO-Behind (B). This extensive study is obtained in 
the light of the recently proposed scenario for the interpretation of the nanodust- 
impact-associated potential pulses measured between the spacecraft main body 
and the boom antennas (Zaslavsky et al., 2012; Pantellini et al., 2012a, b). This 
scenario was enabled by the detailed analysis of the Time Domain Sampler 
(TDS: Zaslavsky et al. 2012) data. 

2. Dust Impact Detection by STEREO/WAVES LFR 
2.1. Description of the Instrument 

The STEREO- A and STEREO-B spacecraft travel close to the Earth's orbit, 
STEREO-A preceding Earth and STEREO-B trailing behind at angular dis- 
tances increasing by 42 degrees per year. The STEREO/WAVES or S/WAVES 
instrument (Bougeret et al., 2008) is the radio and plasma wave experiment 
aboard the STEREO spacecraft. Its sensors consist of three, mutually orthogo- 
nal, six- meter long monopole antennas, which we label X, Y, and Z (Bale et al., 
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2008). The S/WAVES instrument produces power spectra of electric field fluctua- 
tions (Low Frequency Receiver: LFR, and High Frequency Receiver: HFR), along 
with limited samples of the raw voltage timeseries, also called waveforms ( Time 
Domain Sampler: TDS). In this study, we will focus on the LFR data measured 
when high- velocity nanodust particles impact the spacecraft. An extensive study 
of the TDS data during dust events was made by Zaslavsky et at. (2012). In 
dipole mode, LFR measures the difference between the voltages of the X and Y 
antennas about every 40 seconds, which is the typical time resolution for both 
STEREO spacecraft. The LFR frequency range is divided in three frequency 
bands of two octaves: A, B, and C. For each band, the acquisition time of the 
power spectrum is inversely proportional to the frequency at the center of the 
band. Therefore, band A (2.5 kHz - 10 kHz) has the greatest acquisition time 
(406 ms, respectively 4 times and 16 times greater than for bands B and C), 
which makes it the most sensitive to dust impacts, due to their non-stationary 
behavior and to the automatic gain control. 

2.2. Analysis of the Dust Impact Signal 

For dust grains impacting the spacecraft, with rate N, voltage pulse maximum 
amplitude SV, and rise time r (Zaslavsky et ai, 2012, Figure 2), the theoretical 
power spectrum is (Meyer- Vernet, 1985) 

V? « 2(N5V 2 uj- 2 (1 + wV)" 1 ) (1) 

at frequencies [/ = lo/2tt] much greater than the pulses' inverse decay time 
(itself much greater than r), where the angular brackets stand for averaging over 
the pulses detected during the acquisition time. This produces a power spectrum 
varying as /~ 4 at frequencies greater than 1/2ttt, and as f~ 2 at frequencies lower 
than 1/27TT. Comparing Equation (1) to the LFR band A measurements, it is 
possible to determine the rise time [r] and the amplitude [(N5V 2 )]. To do so, it 
is necessary to evaluate the ubiquitous thermal noise (Meyer- Vernet and Perche, 
1989), which decreases less steeply than the spectrum produced by dust impacts. 
The relative intensity of signals from dust and thermal noise is such a strong 
function of frequency that band C, which responds to the highest-frequency sig- 
nals, only detects thermal noise, even during dust impacts (Meyer- Vernet et al., 
2009b). It is noteworthy that contrary to TDS, LFR cannot measure individual 
impacts (Meyer- Vernet et ai, 2009b), due to: i) the Automatic Gain Control 
(AGC), which determines the receiver gain (Bougeret et ai, 2008) and which 
needs multiple dust impacts to record them, ii) the fact that LFR integrates 
all of the impacts occurring during its acquisition time, and iii) the wavelet-like 
transform (Sitruk and Manning, 1995), which is used on board to compute the 
signal power spectral densities. 

Figure 1 shows the power spectrum measured in Band A by STEREO-A and 
the result of our fitting procedure in the case of dust impacts (upper spectrum) 
and when measuring the ubiquitous plasma quasi-thermal noise (lower spec- 
trum). This figure illustrates the difference in power level and spectral index 
between the plasma thermal noise and the detection of nanodusts at these fre- 
quencies. In the case of dust impacts, we fit Equation (1) to the power spectrum 
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Figure 1. Power spectrum measured on Band A by STEREO-A on 10 October 2009 at 5:12 
(upper spectrum) and 7:25 (lower spectrum). The upper spectrum is typical of dust measure- 
ment, whereas the lower one corresponds to the plasma quasi-thermal noise measurement. The 
lines correspond to the best fit, in case of dust (red solid line) and plasma thermal noise (blue 
dashed line). 



measured in Band A to obtain the two free parameters: the rise time r and the 
amplitude (N5V 2 ). 

Previous studies have shown that on STEREO-A the X-antenna is in the best 
geometric configuration to directly detect the plasma cloud from fast nanometer- 
scale dust impacts through the perturbation of the photoelectron population 
around the antenna (Pantellini et al., 2012b; Zaslavsky et al, 2012). Conse- 
quently, the voltage on the X-antenna is typically 20 times larger than on the 
Y antenna (Zaslavsky et al, 2012). The power spectrum measured by LFR in 
Band A (based on the difference in voltage between the X- and Y-antennas) 
is then much higher than the thermal noise, at least two orders of magnitude 
higher. The voltage pulse on the X-antenna is given by (Zaslavsky et al., 2012) 

SV X « ^ (2) 

where T w 0.5 and L = 6 m are the antennas gain and length, respectively (Bale 
et al., 2008), I is the length of antenna within the plasma cloud, and T = 2.5 eV 
is an effective temperature, of the order of the photoelectron temperature as 
implied by the physical process at the origin of SVx (Zaslavsky et al., 2012; 
Pantellini et al, 2012b). 

On the other hand, on the Y- antenna the voltage pulse is 6Vy ~ TQ/C SC , 
where C sc is the spacecraft body capacitance which is about 200 pF (Zaslavsky 
et al, 2012), and the charge Q available in the cloud. The charge [Q] depends 
on mass, speed, angle of incidence, as well as grain and target composition. It 
can be determined semi-empirically with a large uncertainty with (McBride and 
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Table 1. Averaged rise time values year by year and for the whole data 
set for both STEREO spacecraft. N is the number of spectra with de- 
tectable dust impacts, and p is the percentage of the returned data that 
N corresponds to. 







STEREO-A 




STEREO-B 


uc 


N 


p[%] 


t ± At [/is] 


N 


p[%] 


t ± At [/is] 


2007 


158 164 


20.5 


39 ± 15 


29 008 


3.8 


40 ± 22 


2008 


287 705 


36.1 


39 ± 8 


3 030 


0.4 


43 ± 22 


2009 


124 767 


15.6 


40 ± 9 


36 447 


1.5 


39 ± 16 


2010 


79 394 


9.9 


39 ± 12 


1311 


0.2 


44 ± 35 


2011 


1001 


0.1 


42 ± 32 


10 977 


1.4 


39 ± 23 


all 


651 031 


16.4 


39 ± 11 


80 773 


2.0 


39 ± 20 



McDonnell, 1999) 

Q « 0.7m 1 ' 02 1; 3 ' 48 (3) 

where Q is expressed in Coulombs, m in kilograms, and v in kms -1 . 

However, on STEREO-B neither the X-antenna nor the Y-antenna are close 
enough to the preferred impact zone of fast nanometer grains, which is found to 
be close to the Z-antenna (Zaslavsky et at, 2012). Consequently, both antennas 
see a voltage pulse given by 8Vy ps TQ/C sc , so that LFR can observe a dust 
signal only in the rare cases of impacts close to the X-or Y-antennas. Therefore 
on STEREO-B, LFR does not allow accurate measurements of the flux of the 
nanodust, but it still provides information on the rise time [r] in these rare cases. 



3. Nanodust Characteristics at 1 AU 
3.1. Rise- Time Measurements 

Our analysis of the power spectra measured by LFR allows to us determine the 
rise time of the pulses produced by dust-grain impacts. Since LFR does not 
measure individual impacts but integrates over all impacts occurring during the 
acquisition time, the rise-time [r] obtained in our study is an average value. 

Table 1 gives the rise-time values, and Figure 2 presents the histograms of the 
rise-time measured by STEREO-A (panel (a)) and STEREO-B (panel (b)), for 
the whole data set and year by year. Both Figure 2 and Table 1 show that the 
rise time measured by LFR is stable in time and similar for both STEREO-A 
and B, with a value of t = 39 ± 12 yus. 

The rise times measured by LFR come from the integration of all of the 
impacts detected during the integration time. Nevertheless, the individual mea- 
surements of the rise time made using TDS show a similar most probable value 
around 40 /us. The difference between these two measurements is a lack of value 
of t greater than 100 ^s in the LFR data compared to TDS. For STEREO-A, 
the greatest values of r measured by TDS are associated with the largest pulse 
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Figure 2. Histograms of the rise time measured by STEREO-A (panel (a)) and STEREO-B 
(panel (b)), for all data (black solid line) and year by year (red long-dashed line: 2007, green 
dotted-dashed line: 2008, blue dashed line: 2009, cyan dotted line: 2010, magenta solid line: 
2011). 



amplitude [Vjc], corresponding to the largest dust particles. It is not surprising 
that the spectral analysis of LFR is not sensitive to the large dust particles, 
which are the less numerous ones and contribute less to the integral (Equation 
(1)) because of the large r in the denominator. 

The fact that the measured rise times are stable with time and similar for 
both spacecraft is in agreement with the recently proposed scenario for the inter- 
pretation of dust-impact-associated voltage pulses on boom antennas (Zaslavsky 
et al., 2012; Pantellini et al., 2012a, b). In this scenario, the rise time is determined 
by the photoelectron cloud around the antenna, and can be evaluated to be 
proportional to the inverse plasma frequency of this photoelectron cloud. Since 
there are not noticeable changes in the radial distance between the Sun and the 
spacecraft, the plasma frequency of the photoelectron cloud is constant. Then, 
the value of r is a function of the spatially varying photoelectron density [n pn ] 
inside the extended plasma cloud, which is difficult to quantify (Pantellini et al., 
2012a). From our observed value of r, the "mean" value of n p h in the cloud is 
about 10cm~ 3 . This scenario also implies that the expansion time of the plasma 
cloud from the impact zone to its maximum size (i.e. when the plasma cloud 
density reaches the value of the ambient solar-wind density) is greater than the 
measured value of r. 

3.2. Determination of the Flux from LFR Measurements 

The parameters [r and (NSV 2 )] obtained from the fitting of the LFR spectra 
in Band A of STEREO-A enable us to infer the nanodust flux observed at 1 
AU. Considering a dust grain impacting at a distance r from the antenna, and 
a spherical expansion of the impact cloud, the length of the antenna within 
the cloud [I] is given by / = \/Rc ~ 7-2 wncre Rc ~ (3Q/47ren a ) 1 ^ 3 is 
the maximum size of the plasma cloud. The maximum value of the distance 
r for each dust impact is the minimum value between Rq and i?sCi the order 
of magnitude of the spacecraft size (Zaslavsky et al., 2012). Then assuming a 
cylindrical symmetry around the antenna, so that the probability that an impact 
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occurs at a distance between r and r + dr from the antenna is proportional to 
lixrdr, we can write 



(N5V 2 /t 2 ) 



f(m) 



min(ij c ,ij sc ) fiyl 



2itrdrdm 



(4) 



where /(ro)dm = fom 1 dm is the flux of particles of mass between m and 
m + dm. Note that the interplanetary dust models would give /o = 7.4 x 

m max ] (Grun 



10- 19 kg 



; s 1 , and 7 ~ 1 1/6 in the mass range [m r 



et al, 1985; Ceplecha et al, 1998). 

To obtain the nanodust flux from Equation (4), we need the average value 
of the squared impulse voltage produced by a cloud of charge Q: SV 2 . Using 
Equation (2) for 8V and considering both cases of a plasma cloud larger and 
smaller than the spacecraft effective size 



r 2 T 2 R? 



(6V 2 



L 2 



2R s 

(4 - %) 



if Rc < Rsc 
if Rc > Rsc 



(5) 



In the mass range considered in this study and using Equation (3), we have 
R c w Km 1 / 3 , with K w (3 x 0.7w 3 - 5 /47r« a e) 1/3 ~ 5 x 10 6 mkg _1/3 , for typical 
values of the ambient solar- wind density n a , and impact speed v = 300 km s -1 
(Mann et ai, 2010), which is close to the solar wind velocity because of the pick- 
up-ion-like mechanism. For STEREO-A, the effective sensitive area is i?| c « 
0.7m 2 (Zaslavsky et al., 2012). The mass for which i?sc = Rc is m R S c = 
1.6 x 10 _20 kg. This mass is of the same order of magnitude of the maximum 
dust mass detected, m max m 2 x 10~ 20 kg, corresponding to the mass range 
measured by TDS (Zaslavsky et at, 2012). 

Assuming that r is a constant, according to Section 3.1, we integrate Equation 
(4) using Equation (5) 



5V 2 /r 2 ) ~ J nR 2 c ( g 



#2 m - 7 +5/3 



Km 



-7+7/3 1 mR 



2i? 2 c (- 7 + 7/3) 



-7 + 5/3 2(- 7 + l) 




(6) 



Assuming that m mln <C m^ s 
becomes 



?7i. 



and taking 7 = 11/6, Equation (6) 



rr 



(NSV 2 /r 2 ) « fonK'mg, L I ^ ) , 



2 



which only depends on the mass tor sc 
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Figure 3. Variation with time of the cumulative flux of particules of mass greater than 
10 _20 kg measured by STEREO- A/ WAVES LFR between 2007 and 2011. The color scale of 
the cumulative flux distribution has been chosen to emphasize the daily most-probable value 
of the flux. The dotted gray line corresponds to the daily mean of the measured nanodust flux. 

3.3. Nanodust Fluxes at 1 AU 

Figure 3 represents the cumulative flux of particles of mass greater than 10~ 20 kg 
(F 10 -2o), measured by STEREO-A between 2007 and 2011 using Equation (7). 
The cumulative flux is used to allow direct comparison with previous studies 
and interplanetary-dust models, and is given by 



The flux has large fluctuations, as first noted by Meyer- Vernet et al. (2009b) and 
confirmed by Zaslavsky et al. (2012). The white periods on Figure 3 correspond 
to time periods when STEREO- A/ WAVES LFR did not measure dust. This is 
expected to be due, at least in part, to geometrical effects. One clue for the 
geometric origin is inferred by the quasi-continuous flux measured by TDS on 
STEREO-B (Zaslavsky et al., 2012). Another indication of the geometrical origin 
of these fluctuations is the dependence of the flux measured by STEREO-A with 
the latitude in ecliptic coordinates. It is noteworthy than this behavior is not 
seen by STEREO-B (D. Malaspina, private communication, 2012). 

In addition to the temporal variations of the nanodust flux at large time scales, 
the high temporal resolution of LFR enables us to study the fast variations of 
the nanodust flux detected by STEREO-A. Figure 4 shows such variations of the 
cumulative flux, F 10 -2o, during 12 January 2007. In this example, 46% of the 
spectra measured by LFR show nanodust impacts, with long continuous periods 
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00:00 06:00 12:00 18:00 24:00 
12 Jan 2007 Time (UTC) 

Figure 4. Variation during 12 January 2007 of the STEREO- A/WAVES LFR measurements. 
The data for which the thermal noise is dominant, so that the dust flux cannot be measured, 
are plotted in gray, whereas the genuine dust flux measurements are plotted in red. The black 
dotted line shows the average cumulative flux [i^g-20] of nanodusts measured on 12 January 
2007. 

of either dust or thermal-noise measurements, and periods of fast alternation 
between both (as between 19:00 and 21:00). These three behaviors take place at 
different time scales, from a few minutes up to days, and the longest continuous 
measurement of dust in our data set is 11.5 days starting on 13 January 2008. 

Figure 5 shows the histograms of F 10 -2o measured by STEREO-A, when LFR 
measures dusts. The histograms are normalized to their respective maximum. 
One can see in Figure 5 that the distribution of measurement is stable at large 
temporal scales, even if the number of measurements can vary strongly from 
one year to another (from 287 706 in 2008 to only 1001 in 2011, see Table 1). 
These measurements agree with the TDS measurement (Zaslavsky et at, 2012), 
with the previous study of LFR measurements (Meyer- Vernet et ai, 2009b), and 
interplanetary dust distribution models, in which F 10 -2o = 0.04 m~ 2 s _1 (Grun 
et ai, 1985; Ceplecha et at, 1998). The most-probable value of the flux is not 
much larger than the instrumental limit of detection, leading to a very asymmet- 
ric distribution. Even if LFR could measure smaller flux, it is noteworthy that 
flux lower than 2 x 10~ 3 m~ 2 s _1 cannot be measured due to the plasma thermal 
noise, as shown in Figure 4. Since the geometrical effects, like projection effects 
including the flow direction of the dust and the orientation of the spacecraft, 
on dust detection are not fully understood at this time, it is not possible to 
determine the origin of the large difference of dust measurements from one year 
to another. 
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2007, green dottcd-dashcd line: 2008, blue dashed line: 2009, cyan dotted line: 2010, magenta 
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4. Conclusions 

Even if not designed to do so, S/WAVES LFR allows accurate measurements 
of the interplanetary nanodust flux at a high temporal resolution. Nevertheless, 
the detection mechanism needs one of the antenna booms used by the receiver 
to be inside the plasma clouds created by the dust impacts. Consequently, only 
STEREO-A can be used to determine the nanodust flux with LFR. However, 
complementary measurements by TDS allow one to study dust with STEREO- 
B. The complementarity of the time- and frequency-domain measurement of 
nanodust using antenna booms is not limited to this point. Both TDS and LFR 
instruments measure nanodust signals with a similar rise time of 40 /is and a 
similar flux, but with complementary advantages in dust flux measurements. 
TDS measures individual nanodust impacts in monopole mode using the three 
antennas, allowing dust measurement on both STEREO spacecraft but with 
a limited telemetry, whereas LFR gives a continuous survey at a high time 
resolution. 

Geometrical effects on nanodust detection using S /WAVES have been pointed 
out in this study and previous ones, but due to the very complex geometry of 
the spacecraft from the point of view of the antenna booms and in the moving 
frame of the dust, we do not yet fully understand the observed variations of 
the nanodust flux. Nevertheless, this continuous survey of the nanodust in the 
interplanetary medium at high time resolution allows the ongoing study of the 
correlation between the local solar- wind plasma and the properties of the dust. In 
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addition to the TDS measurements, this study should lead to important insights 
on dust detection mechanisms for both STEREO and nanodust acceleration by 
the magnetized solar wind. 

The spectral analysis of more than 700 000 samples shows that the rise time 
of the signal created by dust impacts on the spacecraft is stable, in agreement 
with the recently proposed scenario for the interpretation of the dust impact 
associated voltage pulses on boom antennas (Zaslavsky et al, 2012; Pantellini 
et at, 2012a,b). This enables us to define the spectral and temporal resolution 
and domain required to measure accurately dust with an antenna boom for 
future instruments, such as Fields on Solar Probe Plus. 

From the insight on the detection mechanism, we have updated the previous 
determination of the dust flux measured by LFR (Meyer- Vernet et al., 2009b). 
With our improved interpretation, we found a flux of the same order of mag- 
nitude, but without requiring estimates of unknown parameters, such as the 
expansion velocity of the cloud or the minimal mass of impacting dust. The 
main uncertainty of our measurements comes from the relation between the 
released charge, the mass, and the impact speed (Equation (3)) for nanodust. It 
is noteworthy that this limitation is inherent to all measurements of dust in this 
size range, and that ongoing and future improvements of dust accelerators in 
laboratory experiments should resolve this issue in the future. In this study, we 
have assumed that the power-law index of the interplanetary dust distribution 
models (Grun et al., 1985; Ceplecha et al., 1998) is valid in the nanometer range, 
in other words that the nanodust flux behaves as if these particles were created 
by collision equilibrium. Further analysis of the TDS measurements of nanodust 
could infirm this, leading to an observed distribution of the nanodust flux. In 
such a case, the numerical values of the flux that we give in the present study 
will have to be changed by updating the value of 7 in Equation (6). 
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